Molecular production at a wide Feshbach resonance in Fermi-gas of cooled atoms. 
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The problem of molecular production from degenerate gas of fermions at a wide Feshbach reso- 
nance, in a single-mode approximation, is reduced to the linear Landau-Zener problem for operators. 
The strong interaction leads to significant renormalization of the gap between adiabatic levels. In 
contrast to static problem the close vicinity of exact resonance does not play substantial role. Two 
main physical results of our theory is the high sensitivity of molecular production to the initial 
value of magnetic field and generation of a large BCS condensate distributed over a broad range of 
momenta in inverse process of the molecule dissociation. 



In recent years there have been numerous achievements 
in the area of ultra-cold atomic physics. The major ex- 
perimental tool for it is the use of the Feshbach reso- 
nances (FR) [H, [3, [3, i, H, H, 0|, which occurs when the 
energy of a quasibound molecular state becomes equal 
to the energy of two free alkali atoms. The magnetic- 
field dependence of the resonance allows precise tuning 
of the atom-atom interaction strength in an ultracold gas 
[H. Moreover, time-dependent magnetic fields can be 
used to reversibly convert atom pairs into weakly bound 
molecules i, i, E El E [13, Q [H ■ This technique 
has proved to be extremely effective in converting degjen- 
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erate atomic gases of fermions ^, ^ 
and bosons [3, 113, [l3] into bosonic dimer molecules. 

The Feshbach resonance proceeds at sufficiently strong 
magnetic fields so that electronic spins are polarized. The 
collisions in s-channel of such Fermi-atoms is possible 
only if they have different states of nuclear spins. In a 
typical experiment [1, 17| an admixture of atoms of ^"if 



with the same total atomic spin 9/2 but different spin 
projection quantum numbers —7/2 and —9/2 was used. 

Theoretical works on the molecular production can be 
roughly divided in two categories. The first is a phe- 
nomenology suggesting that the pairs of molecules per- 
form independently the Landau-Zener (LZ) transitions 
6] . Therefore the total number of molecules in the end 
of the process is the LZ transition probability multiplied 
by the number of pairs. The most problematic in this ap- 
proach is what should be accepted for the LZ transition 
matrix element A (further we call it the LZ gap). Direct 
calculation of the transition probability from a micro- 
scopic Hamiltonian up to the 4-th order in the interaction 
constant shows that, in contrast to the assumption of 
phenomenological works, the many-body effects are very 
essential. Another c ateg ory includes the works based on 
a simplified model [2]J, in which molecules have only 
one available state mimicking the condensate [13, l24| . 
Though the numerical works of this category displayed a 
reasonable temperature dependence, they could not give 
a clear physical picture and detailed dependencies on pa- 
rameters of the problem. The series of semi-analytical 
works by Pazy et al. [H HI, IH, [27 1 were based on two 



In this letter we consider the process of molecule pro- 
duction from fermi-gas of atoms after the FR is swept 
across the Fermi sea. The accepted model is valid un- 
der assumption of strong interaction equivalent to the 
condition of wide resonance [S^l . We derive the closed 
equation for such process. We show that the problem 
of molecular production from degenerate gas of fermions 
in a single-mode approximation can be reduced to the 
linear LZ problem for operators. In this respect our re- 
sult agrees with the conjecture of the phenomenological 
theories, but the strong interaction leads to significant 
renormalization of the LZ gap, which occurs to be in- 
dependent on the fermi-gas density. Our results display 
a significant dependence of the molecular production on 
the initial state preparation. At the inverse transforma- 
tion of the BEG molecular gas into atomic gas the latter 
appears in the state with strongly developed BCS con- 
densate directly after magnetic field sweeping. 

Our starting point is the Timmermanns et al. Hamil- 
tonian [3l| 



H = ^(e(t) + €p)al„ap„ + ^ ujpblbp ■ 



(1) 



contradicting assumptions as we show later. 



describing fermionic atoms, created by aj,^ with momen- 
tum p, spin a =t, i which distinguishes between the two 
internal states of the atoms, and kinetic energy , that 
are coupled to diatomic bosonic molecules with kinetic 
energy ujg created by &J. The position and the width 
(molecular lifetime) of the FR are respectively controlled 
by the bare time-dependent detuning energy e{t) and the 
coupling constant g, the former experimentally tunable 
by magnetic field. The Hamiltonian ([1]) neglects nonreso- 
nant atom-atom and molecule-molecule interactions that 
near a FR are subdominant to the resonant scattering. 

In this letter we use the word "atoms" for both, un- 
coupled atoms and those inside molecules, and the word 
"fermions" for atoms that are not bound in molecules. 
Correspondingly, we denote the number of atoms N and 
the number of fermions Np. 

First we demonstrate that the condition of wide res- 



2 



onance allows to neglect kinetic energy of fermions. In- 
deed, the calculations in weak coupling approximation 
(20| show that the characteristic interaction energy is 
g\pn,, where n = N jV is the density of atoms. If g\fn ex- 
ceeds the Fermi energy of fermions = ft^(37r^ri)^/'^/2m 
the fermion dispersion can be neglected. This require- 
ment is equivalent to: 



{2mg/h''f 

(37r2)4/3 



-1/3 » 1 



(2) 



which coincides with the definition of the wide resonance 
[22,]. The same criterion allows to neglect the dispersion 
of molecules. This approximation has been used earlier 
2^ 2^ and numerical proof of its validity has been re- 
ported in [27 1 . 

In what follows we assume that only molecules with 
zero momentum &J^o produced, neglecting molecules 
Vq^o excited above the molecular condensate (the single- 
mode approximation). This assumption was justified for 
the equilibrium state in the case of wide resonance [2^ 
and used in the most of theoretical works on the dynam- 
ics of transition |27| . |29| . We thus replace bq in 
Eq. II]) by b^Sq^ = bSq,o. With this simplification one 
notices that only a pair of fermions with opposite mo- 
menta and spin can be converted into a molecule. The 
two assumptions allow to simplify the Hamiltonian ([1]) 
and to solve the problem exactly. 

Let us introduce spin operators Sz , S- and S+ — Sl_ 
as follows: 



■pi 



-Pi 



1) -S- = X!"pT«-pi; (3) 



with the standard commutation relations [5*+, 5-] = 2Sz, 
[Sz,S+] = S+. Neglecting fermion and Boson kinetic 
energy one can rewrite the Hamiltonian ([1]) in terms of 
spin operators: 



H = 2e{t)Sz 
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(bS, 



(4) 



This Hamiltonian commutes with the operator Q — 
Sz + b'^b equivalent to the total number of atoms N 
and additionally with the square of total spin operator 
— Si + S+S- — Sz- It follows from the first equation 
© that Sz = {Np - Ns)/2, where Np is the number 
of fermions and Ns is the number of available fermionic 
states. Thus, Q ^ {N - Ns)/2. Since Ns:$> N ^ Np we 
can approximate Sz ~ —Ng/2. 

The Heisenberg equations of motion are: 



hb- 



-igS^ 



hS- = -~2ie{t)S- + 2i0Sz (5) 



Generally these equations are non-linear. However, in 
the wide-resonance approximation Sz — —Ns/2, they be- 
come linear. Eliminating S*-, we arrive at an ordinary 
linear differential equation for the operator b: 



tiH + 2ine{t)b + A% ^ 



(6) 



where A = g\fnl and = Ns/V is the density of avail- 
able states. Linear equation ([6]) turns into the parabolic 
cylinder equation if e{t) is a linear function of time. In 
the LZ theory it describes the evolution of the amplitude 
to find the system in one of its two states. The role of the 
LZ gap is played by the value A — g^/ri^, which strongly 
exceeds g\/n. The characteristic time during which the 
LZ transition takes place Tpz is determined by require- 
ment that the instantaneous frequency 2e{TLz) becomes 
equal to A. Thus, the characteristic value of e{t) is equal 
to A = g^/nl ^ ep. 

The model neglecting the dispersion is valid until the 

2 

kinetic energy ^ <^ g^fnl. The value ps limiting avail- 
able states is determined by equation: 



Vs 

2m 



5V"- 



(7) 



and the density of available states is associated with 
the limiting momentum by a standard relation: 



37r2 



2\6 



(3^2)4 



(8) 



Both ris and A are independent on the density of fermions 
n. The condition of the wide resonance ([2]) requires that 
Ns ^ N thus justifying the accepted approximation. In 
a series of works by Pazy, Tikhonenkov et 
the authors neglected dispersion and put N. 



N. 



2i,|23, 

As it 

was demonstrated above the neglecting of dispersion is 
justified only if the condition of the wide resonance ^ is 
satisfied. But the same condition ensures that Ng ^ N. 
Thus, their two fundamental assumptions physically con- 
tradict each other, though their model is mathematically 
consistent. Exact quantum solution of the same problem 
was recently found by Altland and Gurarie [l^l • 

The value g can be extracted from the experimental 
data on the dependence of the scattering length on mag- 
netic field near the PR [IBl using a well-known relation: 
g — h^y 4H(a — ao)e/m. On the other hand it can be 

estimated theoretically as (7 ~ ehf\/a^-, where ehf is the 
hyperfine energy and oq is the radius of the molecule far 
from the resonance. Both these estimates give for 
g 10"28erg x cm^/2 ^nd A = O.OOGfsT. The sweep- 
ing of magnetic field proceeds with amplitude few hun- 
dred Gauss. It corresponds to energy scale about Q.03K, 
larger than A. 

The first important conclusion is that the dynamic LZ 
problem is simpler than the static equilibrium one. In the 
latter problem most interesting and intriguing phenom- 
ena take place at very small e < ep. In this range of en- 
ergy the s-scattering amplitude reaches its unitary limit 
and scattering length changes sign resulting in strong 
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BCS coupling and complex behavior of the atomic and 
molecular densities. It is not the case in the dynamic 
problem since the effective interval for transformations 
of atoms into molecules is A ^ ep. This statement is 
correct for any sweeping rate above relaxation. 

Another important conclusion is that the strong in- 
teraction leads to renormalization of the LZ gap. The 
energy scale which appears in the perturbation theory 
is A^^^ — g^/n In the case of wide resonance 

A = gyfrTs is much larger than A'"-* and docs not de- 
pend on the density of atoms. 

Employing equations ([5]), the general solution of the 
ordinary differential equation reads: 

h{t)^u (t, to) b (to) - igv {t, to) (to) , (9) 
igS- (t) = -ii [t, to) h {to) + i~gv {t, to) S- {to) , (10) 

where u{t,to) and v{t,to) are standard solutions of the 
same equation satisfying initial conditions u{to,to) = 1, 
u{to,to) = and v{to,to) = 0, it {to, to) — 1. These 
solutions have the following properties: 



= A 



y^l = l;u*v + A'^u*v = 
(11) 

The solution (j9|10|) allows to follow the evolution of 
the number of molecules N,n{t) — {b''b){t), the BCS 
condensate amplitude F{t) defined by equation F'^{t) = 
{S+S-){t), and the BCS-BEC coherence factor C{t) = 
{b'^S-){t). In the case when initial coherence factor is 
zero it is given by: 

Nm{t) = \u\^Nm{to)+g^\v\^F^{to) (12) 
~g'F\t) = \u\'NM + ~9'\v\'F\to) (13) 
C{t) = ig-'uu*N„,{to) + ivv*F\to), (14) 



where g = g/VV. Using (fTT|) and summing eqs. (fT2|) 
and (|13p we find that 



NsN^m{t) + F^t) = const, 



(15) 



which is a consequence of the conservation laws. Since 
for any state F'^{t) > 0, if there are no molecules in 
the initial state, their number Nm{t) can not exceed 
a value F'^{to)/Ns at any time. Therefore, the final 
molecule production rate depends on the initial state. 
Below we consider two experimentally most relevant situ- 
ations: no molecules and only molecules and no ferniions 
in the initial state. In both these cases the initial value 
C = {b^S^){tQ) = 0. 

In the case of no molecules in the initial state Nm{ta) = 
general equations p2ll4p simplify to 

Nr,r{t) = g^\v\^F\to): i~gC{t) = -~g''vv* F^{to) 

(16) 

Evolution of F{t) in this case is completely determined 
by Note, that the coherence factor G{t) does not 

remain zero. If at the initial moment there exists both 



condensate of molecules (&) and the BCS condensate 
their time evolution can be obtained from equa- 
tions (|9ll0p by taking averages from both sides. 

The value F'^{to) strongly depends on the initial state. 
We consider first an initial state defined as the filled 
Fermi sphere. The condensate density can be readily 
found for this state: F^ = N/2. Thus, the final av- 
erage number of molecules in this case is very small 
iVm(-t-oo) < N/{2Ns). To produce a reasonable frac- 
tion of molecules it is necessary to have large condensate 
amplitude in the initial state. A natural way to generate 
such initial state is to start with sufficiently small nega- 
tive detuning energy e. The effective dimensionless BCS 
coupling constant reads Xbcs = —vpg^/^, where vp is 
the density of state at the Fermi-energy [28| . This con- 
stant becomes of the order of 1 at | e \^Tep, a value of 
detuning between Fermi-energy and the gap A. When 
detuning becomes less than this value the condensate 
spreads from an exponentially narrow spherical layer near 
the Fermi sphere to a sphere of the radius given by the 
cut-off momentum p,.- The latter is inverse proportional 
to the effective radius of the hyperfine interaction and is 
much larger than ps- Therefore, the molecular produc- 
tion determined by the part of initial BCS condensate 
confined in the sphere of the radius Ps first grows and 
then decreases when the initial detuning decreases. 

Such a strong dependence of the final molecular pro- 
duction on the initial state, in particular on the value of 
the initial magnetic field explains why different experi- 
menters obtain different fractions of molecules in the final 
state even in the adiabatic regime 0, [H, 17, 3]. Note, 
that in experiments in which a significant molecular pro- 
duction was achieved the initial state was indeed close to 
the FR, whereas the final state could be sufficiently far 
from it. 

Thus, in a realistic experimental setup the initial value 
of e is small eo < Fei? <C A and then e increases lin- 
early with time. In this case one can put to = Oi f^nd 
e{t) — et. Equation ([6]) turns into the parabolic cylinder 
equation. Its standard solution u{t, 0) has an asymptotic 
property: |u(cx3,0)p — exp(— 7rA^/2/ie). Employing it 
together with (fTT|) and p^ . we arrive at the following 
number of molecules in the final state: 



Nm{+oo) = F^N;^ [l - exp (-TrAVZ/ie)] 



(17) 



It can be proven that the maximal possible value of F^ 
is NsN/2. It corresponds to complete transformation 
of atoms into molecules in the adiabatic regime e — > 0. 
Equation p?)) looks exactly as the LZ transition prob- 
ability multiplied by an effective number of pairs, but 
in contrast to phenomenological theories p!, 'a] and the 
perturbation theory result ^20] the coefficient in front of 
1/e in exponent does not depend on the initial density 
of the atoms. This theoretical prediction can be checked 
experimentally. 
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Finally, we consider an inverse process with no 
fermions, no BCS condensate and only the molecular con- 
densate in the initial state: (6) (—00) — ■\/iV/2 and sweep- 
ing of the magnetic field in opposite direction. Then in 
the end the condensate density is determined by the LZ 
value: (6)(+oo) = ^iV/2cxp(— 7rA^/2?ie), whereas the 
absolute value of the BCS condensate amplitude {S-) 
can be found from the conservation law (1151): 



N,N 



1 " exp — 



(18) 



This result has a clear physical interpreta- 
tion. It corresponds to [l — exp (— 7rA^//ie)] 
Cooper pairs distributed with identical probability 
w = 2^ [1 ~ 6xp (— 7rA^/?ie)] between Ng available 
states. Then the modulus of the pair amplitude at a 
fixed state is (ap|a_p|) = \/w. If all these amplitudes 
have the same phase, the total condensate amplitude 
is equal to {S-) = Nsy/w, which is equivalent to 
equation (|18p . This result has experimentally verifi- 
able consequence. Indeed, the pair created after the 
sweeping of magnetic field has the size in real space 



'^pair — ^ 



It means that the pair is a compact 
formation well separated from other pairs and, therefore, 
can be rather qualified as a quasimolecule. In contrast to 
real molecules the quasimolecules have parallel electron 
spins. The second their peculiarity is their instability: 
after the sweeping of magnetic field stops, they decay 
during the relaxation time. The latter is rather long 
since the pair fermion collisions do not produce energy 
relaxation. The experimental estimate for the relaxation 
time is in the range of seconds. Therefore, it seems quite 
feasible to switch off the trap for much shorter time 
and observe the correlations of momenta and spins in 
runaway particles. The prediction of our theory is that 
the correlation prefers opposite velocities and parallel 
spins in the range of energy up to A. 

In conclusion, we considered the molecule formation 
and dissociation in cooled Fermi-gas when magnetic field 
is swept across the wide Feshbach resonance. It was 
demonstrated that in this situation the fermion kinetic 
energy is negligible. The resulting molecular produc- 
tion from initial fermions is described by LZ-like for- 
mulae with strongly renormalized LZ gap independent 
on the density of initial fermions. However, the molecu- 
lar production strongly depends on initial value of mag- 
netic field. In the inverse process of molecular dissocia- 
tion immediately after the sweeping stops, there appear 
Cooper pairs with parallel electronic spins and opposite 
momenta homogeneously distributed within a sphere in 
the momentum space, whose radius ps is much larger 
than the Fermi momentum. Another experimentally veri- 
fiable prediction is independence of the coefficient in front 
of 1 /e in LZ exponents entering the equations for molec- 



ular production (|17p and the BCS condensate amplitude 
p8)l on the initial density of atoms (molecules) . 
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